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Abstract. 

We consider in detail the problem of gauge dependence that exists in relativistic 
perturbation theory, going beyond the linear approximation and treating second and 
higher order perturbations. We first derive some mathematical results concerning 
the Taylor expansion of tensor fields under the action of one-parameter families (not 
necessarily groups) of diffeomorphisms. Second, we define gauge invariance to an 
arbitrary order n. Finally, we give a generating formula for the gauge transformation 
to an arbitrary order and explicit rules to second and third order. This formalism 
can be used in any field of applied general relativity, such as cosmological and 
black hole perturbations, as well as in other spacetime theories. As a specific 
example, we consider here second order perturbations in cosmology, assuming a flat 
Robertson- Walker background, giving explicit second order transformations between 
the synchronous and the Poisson (generalized longitudinal) gauges. 
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1. Introduction 

The perturbative approach is a fundamental tool of investigation in general relativity, 
where exact solutions are most often too idealized to properly represent the realm 
of natural phenomena. Unfortunately, it has long been known that the invariance 
of general relativity under diffeomorphisms (two solutions of the Einstein equation 
are physically equivalent if they are diffeomorphic to each other) makes the very 
definition of perturbations gauge dependent jim, ^, ^ ||, ^, |^, ||, |^ . A gauge choice is 
an identification between points of the perturbed and the background spacetimes, and 
generic perturbations are not invariant under a gauge transformation. This "gauge 
problem" has been widely treated in linear theory, but what if one wants to consider 
higher order perturbations? In particular, how do the latter change under a gauge 
transformation? 

Second order treatments have been recently proposed, both in cosmology [^0|, |ll| 
and compact object theory as a way of obtaining more accurate results to be 
compared with present and future observations. For example, in view of the increased 
sensitivity expected from the next generation of detectors, precise computations of 
microwave background anisotropics or gravitational wave production may require 
going beyond the linear regime. In addition, second order perturbations provide a 
reliable measure of the accuracy of the linearized theory (see, e.g., [|2|). In cosmology, 
a second (or higher) order treatment might prove necessary when dealing with scales 
much smaller than the cosmological horizon, where the linear approximation can either 
be not accurate enough or just miss some physical effects. Finally, general relativity 
is an intrinsically non-linear theory, thus it is in principle interesting to look at higher 
order perturbations as a tool for exploring non-linear features. Unfortunately, already 
the second order calculations are almost invariably a computational tour de force, and 
a gauge-invariant treatment is not at hand. As a matter of fact, the computationally 
more convenient gauge does not necessarily coincide with the most interesting one, and 
often different authors work in different gauges, but at present a general formalism 
to deal with second and higher order gauge transformations is not available in the 
literature, although some partial results can be found in ||, [l| and in [|| (and 
references therein). 

The aim of this paper is to fill this gap. Thus, while we shall consider the problem 
of gauge transformations from a general geometrical perspective, from the practical 
point of view our main goal is to derive the effect on a tensor field T of a second 
order gauge transformation. To this end, we shall show that the latter is necessarily 
represented in coordinates by 

= + y (efi),,eri) + e^)) + • ■ • , (1.1) 



where ^(i) and ^(2) £^re two independent vector fields, and that, under (1.1), the first 
and second order perturbations of T transform as 

5f =ST+£^,^^To, (1.2) 
S^f = S^T + 2^5,,, ST + £l^^ To + £^,,^ To . (1.3) 

Equation (1.2) is the usual first order result in terms of the Lie derivative of 
the background tensor field Tq along the vector field ^(1), and the second order 



perturbation 5 T is defined by (3.3) and (3.4). In fact, we derive a general formula. 
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equation (4.6), from which the gauge transformations to an arbitrary order n can be 
deduced, although we wiU give the exphcit expression only for transformations up to 
third order. Furthermore, we show that a tensor field T is gauge-invariant to second 
order if and only if it is gauge-invariant to first order, i.e., 

£^To - (1.4) 

for an arbitrary vector field ^, and, in addition, 

£^6T = . (1.5) 

Actually, we generalize the above result giving a condition for gauge invariance to an 
arbitrary order n. 

The plan of the paper is the following. In the next section we give some 
mathematical results concerning the Taylor expansion of tensor fields under the action 
of one-parameter families (not necessarily groups) of diffeomorphisms. This material 
is necessary for the applications to follow; however, the reader interested only in the 
latter can skip all the proofs, as well as the last two paragraphs of section 2.2; section 



2.4 is useful if one wants to translate our results in the language of coordinates. In 
section 3 we give a general discussion of spacetime perturbations and gauge choices, 
giving in particular a precise definition of the fc-th order perturbation of a general 
tensor field. In section 4 we first define total gauge invariance and gauge invariance 
to order n; then we give the generating formula for a gauge transformation to an 
arbitrary order, and explicit rules for second and third order transformations.! In 
section 5 we consider perturbations of a flat Robertson- Walker model, and apply our 
results to the case of the transformations between two specific gauge choices, i.e., the 
synchronous and the Poisson gauges. Section 6 contains a final discussion. 

In general, we shall work on an m-dimensional Lorentzian manifold of signature 
TO — 2. The abstract mathematical notation is used extensively, but sometimes we 
shall make reference to charts; in this case, coordinate indices ^,1/, . . . take values 
from to TO — 1. Units are such that c = 1. Furthermore, for the sake of simplicity 
we shall always suppose that manifolds, maps, tensor fields, etc., are as smooth as 
necessary, or even analytic (see [EOl for an extension to C fields). 



2. Taylor expansions of tensor fields 

Before considering the specific problem of gauge transformations in relativistic 
perturbation theory we need to establish some general results concerning the Taylor 
expansion of tensor fields, to be used later. For functions on H™, a Taylor expansion 
is essentially a convenient way to express the value of the function at some point in 
terms of its value, and the value of all its derivatives, at another point. Of course, 
this is impossible for a tensor field T on a manifold A4, simply because T{p) and T{q) 
at different points p and q belong to different spaces, and cannot thus be directly 
compared. A Taylor expansion can therefore be written only if one is given a mapping 
between tensors at different points of A4. In this section, we study the case in which 
such a mapping arises from a one-parameter family of diffeomorphisms of A4 , starting 
from the simplest case of a flow (i.e., a one-parameter group of diffeomorphisms) and 
then proceeding to generalize it. Let us first establish some notation. 

t A short account of the material covered in sections 2-4 is presented in |17| . 
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Let tf : ^ M he a, diffeomorphism between two manifolds M and M. For each 

T^ip-sM between the tangent 



V(p) 



T*M between the 



p (z AA, Lp defines naturally the linear map Lp^ 
spaces, called push-forward, and the linear map 
cotangent spaces, called the pull-back.^ Using (p~^, we can define also a push-forward 
on T^^ >A/', and a pull-back of T^^p^A/" on TpAi, so that </?*|p and '^*\^(^p'^ turn 
out to be well-clefined for tensors of arbitrary type. We can also construct maps 
and if* for tensor fields, simply requiring that, Vp G A4, 



and 



Hereafter, we drop the suffixes in and ip* 
confusion. 



(2.1) 



(2.2) 

, since there is no real danger of 



2.1. Flows 

Let be a differentiable manifold, and let f be a vector field on Ai, generating a 
flow : IR X — > A4, where (l>{0,p) — p, Vp G A^.f For any given A G H, we shall 
write, following the common usage, 0a (p) := 4>{^,p)i "ip ^ M. Let T be a tensor field 
on M.. The map t/)^ defines a new field 4)\T on A/(, the pull-back of T, which is thus 
a function of A. 

Lemma 1: The field admits the following expansion around A = 0: 



+ 00 



Proof: By analyticity we have 

A^ d^ 
"fc[ dAfc 



= V — — 



T 



where, here and in the following. 



dA* 



(...) := 



dA'' 



(2.3) 



(2.4) 



(2.5) 



and the first derivative is, by definition, just the Lie derivative of T with respect to ^: 



d_ 
dA 



T = lim i {(j)lT - T) £(^T . 



A^O A 

In order to prove ( |2.3| ), it is then sufficient to show that, Vfc, 
k ^' 



dA'^ 



(2.6) 



(2.7) 



I Here, we are following the most common notation, although some authors (see, e.g., [^]|) denote 
the push-forward and the pull-back exactly in the opposite way. 

t In order not to burden the discussion unnecessarily, we suppose that <)> defines global 
transformations of J\A p^ . 



This can be established by induction over k. Suppose that (2.7) is true for some k. 
Then 



hm - 



dA'=^ 



d^ 

dA^ 



T = 



^*-^-dA^- 



T = 



(2.8) 



where t := A — e, and we have used the property that 4i\ forms a one-parameter group: 



0T o 4>e 



□ 



It is worth noticing that (2.3) can also be written [ p3| in the symbolic form <j)\ = 
exp(A£{). 

Equation can be applied to the special case in which the tensor T is just one 
of the coordinate functions on M., x'^. We have then, since (f>*^x^{p) = x^{(f>x{p)), the 
usual action of an "infinitesimal point transformation" extended to second order in A: 



A^ 



(2.9) 



where we have denoted x^{p) simply by x^, and x^((f>\(j>)) by i'' 



2.2. Knight diffeomorphisms 

Let us now suppose that there are two vector fields ^(i) and ^(2) on -M. Separately, 
they generate the flows 4>^^^ and (/j^^^ , respectively. We can combine and to 
define a new one-parameter family of diffeomorphisms ^ -.TRxAi — > A4, whose action 
is given by := 4''xy2 ° ■ Thus, ^'a displaces a point of a parameter interval 
A along the integral curve of ^(ij , and then an interval A'^/2 along the integral curve of 
?(2) (see figure |l|). For this reason, we shall call it, with a chess- inspired terminology, 
a knight diffeomorphism. 

This concept can be immediately generalized to the case in which n vector fields 
• ■ ■ J C(n) are defined on Ai, corresponding to the flows 4>^^\ . . . , Then we 

define a one-parameter family ^ : IR x ^ of knight diffeomorphisms of rank n 

by 



(2.10) 



and the vector fields ^(i), . . . will be called the generators of ^. 

Of course, ° ^\ 7^ '^a+x', consequently, Lemma 1 cannot be applied if we want 
to expand in A the pull-back ^f^T of a tensor field T defined on Ai. However, the 
result can be easily generalized. 



Lemma 2: The pull-back ^/'^T of a tensor field T by a one-parameter family of knight 
diffeomorphisms \1/ with generators ^(i) , . . . , 'C(fe) i ■ ■ ■ can be expanded around A = as 
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follows: 



+00+00 +00 \h+2l2 + ---+kh + --- 

^'*T = y y . . . y . . . - — — - — — — — - — i^'^ ...£''= . . . r . 



(2.11) 



Proof: 



;i=o 

\ ^ \ ^ ^ nh ph Plk rp 



+00 -\-oo -\-oo 



(2.12) 



where we have repeatedly used (O). □ 
The explicit form of (2.11) up to the third order in A is 

^IT = T + A£e^,,T + ^ (4^^ + T 

+ ^ (4i,+3^«a,A(.)+^«(3,)^^+-- - ■ (2-13) 

Equations (2.11) and (2.13) apply to a one-parameter family of knight 
diffeomorphisms of arbitrarily high rank, and can be specialized to the particular 
case of rank n simply by setting ^(j.) = 0, Vfc > n. Their meaning is particularly clear 
in a chart. Denoting, as before, by the coordinates of a point p, and by those 
of 'I'a(p), we have, to order A^, 

= + Ae(% + y + (^2)) +■■■ ■ (2-14) 

Equation (2.14) is represented pictorially in figure ^ 

Since o'iix ^ *cr+A, and 5*^^ 7^ ^'-a, one may reasonably doubt that forms 
a group, except under very special conditions. This is confirmed by the following 

Theorem 1: The only cases in which '5 forms a group are those for which 
^(fc) = ak^(i), Vfc > 2, with ak arbitrary numerical coefficients. Then, under the 
reparametrization A :— /(A), with /(A) := A + J2k^ auX^ /k\, \l/ reduces to a flow in 
the canonical form. 

Proof: Let us first show that ^(fe) ~ a/cC(i) ; V A: > 2, is a sufficient condition for ^E* to be 

a group. Since Lemma 1 implies (ji'a^ — (f>a}ij, we have 5* a = • • • o 4'al\k/k{ ° ' ' ' ° '/'i^'' = 

Thus, (i) 1-^ o ^-A = o 0^^^ = = 0^ ^ ^^^Yi T = f-^{a + A), and 

(ii) = 0^^)-^ = = with p = /-i(-A). 

To prove the reverse implication, let us consider first the case of a knight 
diffeomorphism of rank two. If ^' has to form a group, then for any A, cr G IR, there 
must exist a r G H such that o ^a ~ ^r, i-e., 

0/,oc^i^K^%^oc^['^^0^^ocl>i^K (2.15) 
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Applying Lemma 1 to ( 2.15| ) we get, to second order in the parameters and for an 
arbitrary tensor T, 

(r - A - a) £^^,^T+ 1 (r^ + - 2At - a^) ^^^T 

+ \{t^ -\^-cj^)£^,,,T+---^Q. (2.16) 

In the limit A, cr ^ 0, one gets from (2.16) that t = A + cr — 0:2 Act — /32A^ ~ 720'^ + • • ■ > 
where a2,/^2,72 are unspecified numerical coefhcients. Substituting back into (2.16) 
we have that 

Act {£^t2) ~ "2£c(i)) T - ^2A^ i^jdjT - 72CT^£^(j,r H = . (2.17) 

It is clear that ( |2.17l ) can be satisfied only if /32 = 72 = and ^(2) — 0^2^(1) • Similarly, 
considering higher rank knight diffcomorphisms, one can show that ^(fe) — Oik£,(i) , Vfc.D 

The failure of 5" to form a group is also related to the following circumstance. For 
any p e A^, one can define a curve Up : IR ^ 7W by Up(A) := ^'a(p)- However, these 
curves do not form a congruence on Ai. For the point Up{X), say, belongs to the image 
of the curve Up, but also to the one of u^^^y^ , which differs from Up when at least one of 
the ^(fc) is not collinear with since 'u„^(a)(o-) = *£tO*a(p) 7^ *a+<t(p) = ^^(A + ct). 
Thus, the fundamental property of a congruence, that each point of Ai lies on the 
image of one, and only one, curve, is violated. 

Finally, it is perhaps worth pointing out that the result of Lemma 2 cannot be 
written, even formally, as 

fc=0 



with 



because this expression fails to agree with (2.11) for fc > 3. One might try to define a 



vector 77(A) for which an analog of (2.18) holds, but this does not seem very useful or 
illuminating. 



2.3. General case 

Knight diffcomorphisms are of a very peculiar form, and the previous results seem 
therefore of limited applicability. This is, however, not the case, because any one- 
parameter family of diffcomorphisms can always be regarded as a one-parameter fam- 
ily of knight diffcomorphisms — of infinite rank, in general — as shown by the following 

Theorem 2: Let ^:IRxA^^A^bea one-parameter family of diffcomorphisms. 
Then 3 . . . , . . . , one-parameter groups of diffcomorphisms of A^, such that 

*A = ---o0^'l^ o...o0(^l) O0« . (2.20) 
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Proof: Consider the action of ^'a on a function f : Ai ^ M,. A Taylor expansion of 
gives 



fe=0 



The differential operator defined by 
d 



dA 



(2.21) 



(2.22) 



is clearly a derivative, so we can define a vector ^(i) through £^^-^^f ■— ^{i)f ■ Similarly, 



d^ 

dA2 



nf-£lJ 



(2.23) 



is also a derivative, as one can easily check. Thus, we define the vector ^(2) such that 
£i(2)f ■— '^(2)/; so on at higher orders. Hence, we recover (2.11) for an arbitrary 



/. But if (p and ip are two diffeomorphisms on A4 such that Lp* f = ijj* f for every / 
it follows that (p = tp, as it is easy to see. Thus, we establish (2.20). 



□ 



It must be noticed that, although we have supposed so far that maps and fields 
are analytic, it is possible to give versions of Lemmas 1 and 2, and of Theorem 2, that 
hold in the case of ob jects [po| . The main change is the substitution of Taylor 
series like the one in ( 2^ ) by a finite sum of n — 1 terms plus a remainder |Q . The 
meaning of Theorem 2 is then that any one-parameter family of diffeomorphisms can 
be approximated by a family of knight diffeomorphisms of suitable rank. 



2.4- Interpreting the literature: What is what 

The abstract mathematical notation that we have used so far is the most appropriate 
one for the study of gauge transformations in perturbation theory from a general 
point of view. In order to make explicit calculations in special cases of physical 
interest, however, one must introduce a chart. Most of the literature on the subject, 
therefore, is written in the language of coordinates. For this reason, we indicate 
here how to establish a correspondence between the two formalisms. For the sake 
of simplicity in the notation, we shall restrict ourselves to consider the action of the 
pull-back on a vector field, the extension to one-forms and tensors of higher rank being 
straightforward. 

Let therefore {U, x) be a chart of M, with W C an open set and x : U H"* 
given by a; : p I— > {x°{p), x^{p), . . . , a;™~^(p)), Vp £ U. Since the function x'^ : M ^ JR 
is differentiable, we can define the linear map , that associates to a vector on A4 its 
fj,-th component in the coordinate basis defined by the chart (W,x).f 

Consider now a vector field Z and its pull-back Z := '^*^Z, which for every value 
of A is a new vector field on A4. For each point p eU, the components of Z{p) in the 
chart {14, x) are 

Z'^ixip)) - (xi:Z^ {x{p)) , (2.24) 

t Notice that, in tlic language of differential forms, xf^ = dx^*. 
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which becomes, using ( p.2[ ) and the definition of Z, 

Z'^ix{p))^x^,[z{p))=x^iinZ)ip))^{x^o^>-') jZi^xip))) ■ (2.25) 

Now, let us define a new chart {^\{lJ),y)^ with := x^ o 'i'^^. In this way, the 
y-coordinates of the point q := '^\{p) coincide with the a;-coordinates of the point p 
from which q has come under the action of the diffeomorphism: y^{q) — x^{p). We 
have then 

Z^{<P)) = V'i {Z (*a(p))) = ivtZ) {y {^,{p))) = {y>:Z) (y{q)) . (2.26) 

Denoting by Z^^ and Z'^^ the components of Z in the charts {U, x) and ('^\{U), y), 
respectively, we can then write 



Z^^ixip)) - Z'^{y{q)) 



^y^' 
dx^ 



Z-{x{q)) . (2.27) 

x{q) 



That is, the pull-back of Z is characterized by having, in the chart {U, x) at point 
p, the same components that the original vector field has in the chart {'^\{l4),y) at 
point q = ^'a(p). Or, since y{q) = x{p), one can just write Z^{x) = Z'^(x), where 
X simply stands for a point of IR™. This property can also be used to define Z, 
and corresponds to a passive interpretation of the map x, regarded as generating 
a change in the chart on A4 rather than a transformation of A4 {active view). Of 
course, the two viewpoints are equivalent plj l, but the active interpretation is much 
less confusing. 

We end this section by writing down the explicit expression, up to second order in 
A, for the coordinate transformation a; — > y associated to a family of diffeomorphisms 



with generators C(i)jC(2)j ■ • • ■ Calling q :— ^'a(p), we have from (2.14), 

x'^iq) = x^ip) + X^^,^{x{p)) + y (efi)_^(a:(p)) e(i)(x(p)) +er2)(2^(p))) + ■ • • • (2.28) 
By definition, we have also 

y^^iq) ■.= x'^{p)=x'^{q)-X^^,^ix{p)) 

Expanding the various quantities on the right hand side around g, (2.29) becomes 
finally 

y^(<z)=x^(g)-Ae['i)(a:(g)) + y {^^^^ Jx{q)) C^iMq)) ~ ^^,^{x{q))) + ■ ■ ■ . (2.30) 



Equations ( p!28| ) and ( ^3o| ) express the relationship, in the language of coordinates 



between the active and the passive views. Whereas ( 2.28| ) provides us with the 



coordinates, in the same chart {lA, x), of the different points p and q — 'I'a(p), equation 



(2.30) gives the transformation law between the coordinates of the same point q in the 



two different charts (U,x) and {'^\{lA)^y). An equivalent form of the transformation 



(2.30) was already used by Taub in studying the gauge dependence of an approximate 
stress energy tensor for gravitational fields jl^ 



Using (2.30) for the actual computation of the coordinate transformation in (2.27) 



and expanding every term again at second order around x{p), one can derive the 
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components in the chart {U, x) of the puU-back '^*^Z of Z, given in terms of a second 
order expansion formula involving Z and its partial derivatives along ^(i) and ^(2). 
Then, properly collecting the various terms, one can check that this leads to the 
components of the right hand side of (2.13), i.e., to the components of a Taylor 
expansion of ^f'^Z in terms of the Lie derivatives along ^(ij and ^(2) of Z. 

3. Perturbations of spacetime and gauge choices 

In relativistic perturbation theory one tries to find approximate solutions of the 
Einstein equation, regarding them as "small" deviations from some known exact 
solution — the so-called background. The perturbation AT in any relevant quantity, 
say represented by a tensor field T, is defined as the difference between the value T has 
in the physical spacetime, and the background value Tq. However, it is a basic fact of 
differential geometry that, in order to make the comparison of tensors meaningful at 
all, one has to consider them at the same point. Since T and Tq are defined in different 
spacetimes, they can thus be compared only after a prescription for identifying points 
of these spacetimes is given. A gauge choice is precisely this, i.e., a map between the 
background and the physical spacetime. Mathematically, any diffeomorphism between 
the two spacetimes provides one such prescription. A change of this diffeomorphism 
is then a gauge transformation, and the freedom one has in choosing it corresponds to 
the arbitrariness in the value of the perturbation of T at any given spacetime point, 
unless T is gauge-invariant. This is the essence of the "gauge problem," which has 
been discussed in depth in many papers |l|, ||, H, ||, ||] and review articles |, |,|9|. 

In order to discuss higher order perturbations and gauge transformations, and 
to define gauge invariance, we must formalize the previous ideas, giving a precise 
description of what perturbations and gauge choices are. Here we shall mainly follow 
the approach used in references ^, 0, |2l| (cf . also [|4[ |6| ) . 

Let us thus consider a family of spacetime models {{A4, gx^Tx)}, where the metric 
g\ and the matter fields (here collectively referred to as tx) satisfy the field equation 

£[gx,Tx]^0, (3.1) 

and A G M. We shall assume that gx and tx depend smoothly on the parameter 
A, so that A itself is a measure of the amount by which a specific {A4, gx,Tx) differs 
from the idealized background solution {A4,go, tq), which is supposed to be known. In 
some applications, A is a dimensionless parameter naturally arising from the physical 
problem one is dealing with. In this case one expects the perturbative solution to 
accurately approximate the exact one for reasonably small A (see, e.g., [|l2|). In other 
problems, A can be introduced as a purely formal parameter, and in the end, for 
convenience, one can thus choose A = 1 for the physical spacetime, as we shall do in 
section 5. 

This situation is most naturally described by introducing an (to -|- l)-dimensional 
manifold A/", foliated by submanifolds diffeomorphic to A^, so that Af — Ai x M.. 
We shall label each copy of M by the corresponding value of the parameter A. The 
manifold TV has a natural differentiable structure which is the direct product of those 
of Ai and H. We can then choose charts in which x^^ (/i = 0, 1, . . . ,to — 1) are 
coordinates on each leave A4a, and = A. 

Now, if a tensor field Tx is given on each AAx, we have that a tensor field T is 
automatically defined on Af by the relation T{p,X) :— Tx{p), with p e A^a-I In 

t It is worth noticing that tensor fields on Af constructed in this way arc "transverse," in the sense 
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particular, on each A4\ one has a metric g\ and a set of matter fields tx, satisfying 



the field equation (3.1); correspondingly, the fields g and r are defined on Af. 

We want now to define the perturbation in any tensor T, therefore we must 
find a way to compare T\ with Tq. As already said, this requires a prescription for 
identifying points of A^a with those of A4o- This is easily accomplished by assigning 
a diffeomorphism ipx : TV — * A/" such that fxl^^ ■ Mo M-x- Clearly, (px can be 
regarded as the member of a flow tp on Af, corresponding to the value A of the group 
parameter. Therefore, we could equally well give the vector field X that generates 
ip. In the chart introduced above, X™ = 1 but, except for this condition, X remains 
arbitrary. With a slight abuse of terminology, we shall sometimes refer also to such a 
vector field as a gauge. 

The perturbation can now be defined simply as 



ATx := plT\ 



Tn 



(3.2) 



The first term on the right hand side of (3.2) can be Taylor-expanded to get 



where 



ATx 



A 



k=l 



dAfe 



(3.3) 



(3.4) 



X=Q,Mq 



Equation (3.4) defines then the fc-th order perturbation of T. Notice that ATa and 5^T 
are defined on A^o; this formalizes the statement one commonly finds in the literature, 
that "perturbations are fields living in the background." It is important to appreciate 
that the parameter A labeling the various spacetime models serves also to perform the 
expansion ( [3.3| ), and determines therefore what one means by "perturbations of the 
A:-th order." However, as we have already pointed out, there are applications where 
A is, to a large extent, arbitrary. In these cases, the split of ATa into perturbations 
of first order, second order, and so on, has no absolute meaning, because a change of 
A, i.e., a reparametrization of the family of spacetimes, will mix them up. What is 
invariantly defined, is only the quantity ATa, whereas the various 5^T are meaningful 
only once a choice of the parameter has been made. 



Now, we are interested in those cases in which (3.1) is too difficult to solve exactly, 
so that one looks for approximate solutions, to some order n. In fact, we can now 
obtain much simpler linear equations from ( ^.l| ). At first order, differentiating (3.1) 
with respect to A and setting A equal to zero, one obtains pH a linear equation for &g 
and St. At second order, a second derivative with respect to A of (3.1) at A = gives 
an equation of the type 



L[5^g,5^T]^S[5g,5T] , 



(3.5) 



which is linear in the second order perturbations S^g and (5^r, and where the first 
order perturbations 5g, 6t now appear as known source terms. This can obviously be 
extended to higher orders, giving an iterative procedure to calculate A^a and Ata — 
hence gx and t\ — to the required accuracy. 



that their m-th components in the charts we have defined vanish identically. 
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4. Gauge invariance and gauge transformations 

Let us now suppose that two vector fields X and Y are defined on M , such that they 
have X"^ = Y"^ = 1 everywhere. Correspondingly, their integral curves define two 
fiows (p and tp on TV, that connect any two leaves of the foliation. Thus X and Y 
are everywhere transverse to the A4\, and points lying on the same integral curve of 
either of the two are to be regarded as the same point within the respective gauge: ip 
and ip are both point identification maps, i.e., two different gauge choices. 

The fields X and Y can both be used to pull back a generic tensor field T, and to 
construct therefore two other tensor fields f^T and ip^T, for any given value of A. In 
particular, on Mq we now have three tensor fields, i.e., Tq, and 

T^:=^IT\,, T^:=rxT\o. (4.1) 

where, for the sake of simplicity, we have denoted the restriction to TWo of a tensor 
field defined over TV simply by the suffix 0. 

Since X and Y represent gauge choices for mapping a perturbed manifold A4x 
into the unperturbed one Aio, and are the representations, in T W q, o f the 
perturbed tensor according to the two gauges. We can write, using ( |3.2| )-(3.4) and 
Lemma 1, 

+ 00 , J, +00 .J, 

A:=0 ■ fe=0 

+ 00 .J. +00 .J, 

^a'' = E ^'^^ = E fc! ^yT\o = To + A^T, . (4.3) 
fc=o ■ fc=0 



4-.1. Gauge invariance 

If — , for any pair of gauges X and Y, we say that T is totally gauge-invariant . 
This is a very strong condition, because then (4.2) and (4.3) imply that S^T-^ = 5''T^ , 
for all gauges X and Y and for any k. In any practical case one is however interested in 
perturbations to a fixed order n; it is thus convenient to weaken the definition above, 
saying that T is gauge-invariant to order n iff S'^T^ — S^T^ for any two gauges X 
and Y, and Vfc < n. We have then the following {S°T := Tq, ST S^T) 

Proposition 1: A tensor field T is gauge-invariant to order n > 1 iff £^S^T = 0, for 
any vector field ^ on and Vfc < n. 

Proof: Let us first show that the statement is true for n = 1. In fact, if ST-^ = ST^ , 
we have £x-yT\o = 0. But since X and Y define arbitrary gauges, it follows that 
X—Y is an arbitrary vector field ^ with ^"^ — 0, i.e., tangent to M. Let us now suppose 
that the statement is true for some n. Then, if one has also 6"^^T'^\q — S"^^T^\o, it 
follows that £ x-yS^T-^ = 0, and we establish the result by induction over n. □ 

As a consequence, T is gauge-invariant to order n iff Tq and all its perturbations 
of order lower than n are, in any gauge, either vanishing, or constant scalars, or a 
combination of Kronecker deltas with constant coefficients. Thus, this generalizes to 
an arbitrary order n the results of references ||l|, ^, |^. Further, it then follows that 
T is totally gauge-invariant iff it is a combination of Kronecker deltas with coefficients 
depending only on A. 
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^.2. Gauge transformations 

If a tensor T is not gauge- invariant, it is important to know how its representation on 
A^o changes under a gauge transformation. To this purpose, it is useful to define, for 
each value of A S K, the diffeomorphism <i>A : A^o A^o given by 

*A := V-\ ° i>\ ■ (4.4) 

The action of is illustrated in figure We must stress that <i> : IR x M.q so 
defined, is not a one-parameter group of diffeomorphisms on Mq. In fact, ^ 'I'^^, 
and $A+CT 7^ 'i'cr o *&A, essentially because the fields X and Y have, in general, a 
non vanishing commutator, as depicted in figure |4| However, Theorem 2 guarantees 
that, to order n in A, the one-parameter family of diffeomorphisms $ can always be 
approximated by a one-parameter family of knight diffeomorphisms of rank nf (see 
figure ^ for the action of <I>a to second order). 

It is very easy to see that the tensor fields and defined by the gauges <f 
and ip are connected by the linear map <I>^: 

Tl = rxT\^ - {rx^U^\T)l = ^1{^IT\ = ^IT^ . (4.5) 

Thus, Theorem 2 allows us to use (2.11) as a generating formula for a gauge 
transformation to an arbitrary order n: 

^ \^ \ ^ Ph Ph Plk rpX 

••• ••• ...{km ^«(i)^«(2) ■■■^?(.) • 



(4.6) 



To third order, we have explicitly 

A^ 



+ ^ (^|<,, + 3 ^e<.) + ^e(3, )t^ + --- ^ (4-7) 

where ^(i) and ^(2) are now the first two generators of <I>a, or of the gauge 
transformation, if one prefers. 

We can now relate the perturbations in the two gauges. To the lowest orders, this 
is easy to do explicitly: 

Proposition 2: Given a tensor field T, the relations between the first, second, and 
third order perturbations of T in two different gauges are|: 

5r^-5T^ = i:5^^,To; (4.8) 
S^T^ -S^T""^ (i:^,,, -t- £l^^ ) To + 2£^,,, ; (4.9) 

S^T^ - S^T^ = (^£^(3) + 3£^(^)£4,2) -I- i^lij To 

+ 3 + 4,,) ST + 3£5,,/2yX _ (4 

t This result confirms a claim in ]l4| . 

I A second order gauge transformation equivalent to (4.9) has recently been given in psj, sec their 
section III C. 
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Proof: Substitute (4.2) and (4.3) into (4.7). □ 



This result is consistent with Proposition 1, of course. Equation ( |4.8| ) implies that 
T\ is gauge-invariant to the first order iff <£{To — 0, for any vector field ^ on Ai. 



In particular, one must have £5^2)^0 — 0, and therefore (4.9) leads to £^ST — 0. 
Similarly, one has then £^S^T = from (4.10), and so on recursively at higher orders. 

It is also possible to find the explicit expressions, in terms of X and Y, for the 
generators ^(j.) of a gauge transformation: 

Proposition 3: The first three generators of the one-parameter family of 
diffeomorphisms $ are: 

(4.11) 

^i2)^[X,Y]; (4.12) 

C(3) - [2X-y,[X,r]] . (4.13) 

Proof: Substituting (4.2) and (4.3) into (4.7), using the fact that £^^^.^X — 0, and 
identifying terms of first order in A, wc find 

£^^,^To = £y-xT\, . (4.14) 



Since — X™ = and T is arbitrary, we have (4.11). Substituting back, and 
identifying terms of order A^, we have now, similarly, 

£i,.,To = £[x,Y]Tl . (4.15) 



But [X, y]™ — 0, so we obtain (4.12). Analogously, one finds 



^^jgjTo = £[2X-Y,[x,Y]]T\f^ . (4.16) 



Since [2X - Y, [X, F]]™ ^ 0, we get (4.13). □ 



5. An example from cosmology 

As an example of the applications of the gauge transformation obtained, we now 
show how the perturbations on a spatially flat Robertson- Walker background in two 
different gauges are related, up to second order. We shall first consider the metric 
perturbations, then those in the energy density and 4-velocity of matter. Thus in this 
section we choose m = 4, so that the Greek indices fi,!/, . . . take values from to 3, 
and the Latin ones . . . from 1 to 3. 

The components of a perturbed spatially flat Robert son- Walker metric can be 
written as 

5oo = -«(r)2(^l + 22^^^'-)j , (5.1) 



9oi 



+00 -. 



(5.2) 
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(+00 \ +00 

r—l / r—1 



(5.3) 



wheref xi^^^ = 0, and r is the conformal time. The functions 'ip^^\ ^i^K 4>^^\ and xfj^ 
represent the r-th order perturbation of the metric. 

It is standard to use a non-local splitting of perturbations into the so-called scalar, 
vector and tensor parts, where scalar (or longitudinal) parts are those related to 
a scalar potential, vector parts are those related to transverse (divergence-free, or 
solenoidal) vector fields, and tensor parts to transverse trace-free tensors. Such a 
splitting generalizes the Helmholtz theorem of standard vector calculus (see, e.g., 
Q), and can be performed on any spacetime (see, e.g., and references therein) 
imposing suitable boundary conditions. In our case, the shift uj'^p can be decomposed 
as 

c^W^a.c.WII-l-..^'-'^, (5.4) 

where is a solenoidal vector, i.e., o'Wj = 0. Similarly, the traceless part of 

the spatial metric can be decomposed at any order as 

X? - %X^'^)" + S^x^^ + a.x^^ + xl?^ , (5.5) 

where x^'^'"' is a suitable function, x!*^^^ is a solenoidal vector field, and 9'Xij''^ = 0; 
hereafter. 



- -Jy V^ . (5.6) 



Now, consider the energy density /i, or any other scalar that depends only on r at 
zero order: this can be written as 



+00 ^ 



r=l 



For the 4-velocity it^ of matter we can write 

In addition, is subject to the normalization condition u'^u^ = — 1; therefore at any 
order the time component v'^^.^ is related to the lapse perturbation, V'(r)- For the first 
and second order perturbations we obtain, in any gauge: 

- - ^(1) ; (5-9) 

^(2) = - ^(2) + SV-^i) + + «f . (5.10) 

The velocity perturbation u^^^ can also be split into a scalar and vector (solenoidal) 
part: 

t Indices are raised and lowered using S^^ and Sij, respectively. 
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As we have seen in the last section, the gauge transformation is determined by the 
vectors ^(r)- Sphtting their time and space parts, one can write 

^ir) = "''■^ , (5.12) 

and 

= d'p(-) + S-)' , (5.13) 

with did''^' = 0. 
5.1. First order 

We begin by reviewing briefly some weU-known results about first order gauge 
transformations, as we shall need them in the following. From now on, we will drop 
the suffixes X and Y used previously to denote the "old" and "new" gauge choices, 
simply using a tilde to denote quantities in the new gauge. 

For the sake of completeness, we recall here the basic coordinate expressions of the 
Lie derivative along a vector field ^. For a scalar /, a vector Z and a covariant tensor 
T of rank two, these are, respectively: 

£d = fA"" ; (5.14) 

£^z^ = ztr - ; (5.15) 

£iT^. = T^.A" + + C V ■ (5.16) 

Expressions for any other tensor can easily be derived from these. 



5.1.1. General transformation From (4.S), it follows that the first order 
perturbations of the metric transform as 

= 5g^iv + »f«(i)3^u] , (5.17) 

where g^^v is the background metric. Therefore, using (5.16), we obtain the following 
transformations for the first order quantities appearing in (^)-(|5^: 

a' 

V'(i) = '0(1) + "(1) + — «(!) ; (5.18) 

^^^Jf'-<^+^' + <^"^ (5.19) 



^ V2/3(i) --«(!) ; (5.20) 



xg)=x!?+2D.,/3«+rfg+4^; (5.21) 
where a prime denotes the derivative with respect to r. 



For a scalar /i, from (4.8), (5.7), and (5.14) we have 



6fi = 5fj, + /i(o)"(i) ■ (5.22) 
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For the 4- velocity u'^, we have from ( |4.8| ) 

+ . (5.23) 



Usmg (5.15) and (5.8) this gives: 

a' 

= ^"d - - "(1) ' (5-24) 

^(1) = ^(1) " ^(1) " ^(1) ■ (^-2^) 
The 4- velocity is however subject to the constraint (5.9), therefore (5.24) reduces 



to (5.18) 



5.1.2. Transforming from the synchronous to the Poisson gauge Let us now consider 
the particular case of the transformation from the synchronous to the Poisson gauge. 
The synchronous gauge has been the one most frequently used in cosmological 
perturbation theory; it is defined by the conditions goo = ^'^{'''Y' goi = 
p8[ . In this way the four degrees of freedom associated with the coordinate (or 
diffcomorphism) invariance of the theory are fixed. The Poisson gauge, recently 
discussed by Bertschinger jl9|, is instead defined by uj/^'^'^ — Xij = 0. Then, 
one scalar degree of freedom is eliminated from goi (a;*^''-'" — 0), and one scalar and 
two vector degrees of freedom from gij (x^''^" = Xi"^^ = 0)- This gauge generahzes 
the well-known longitudinal gauge to include vector and tensor modes. This gauge, 

(r) (r) 

in which lj,. ' — xlj = 0, has been widely used in the literature to investigate the 
evolution of scalar perturbations Since the vector and tensor modes are set to 
zero by hand, the longitudinal gauge cannot be used to study perturbations beyond 
the linear regime, because in the nonlinear case the scalar, vector, and tensor modes 
are dynamically coupled. f 



Given the perturbation of the metric in one gauge, it is easy to obtain, from ( 5.18 )- 

( F21I) , the gauge transformation to the other one, hence the perturbations in the new 



gauge. In the particular case of the synchronous and Poisson gauges, we have: 

^ + ; (5.26) 



a 



(5.27) 



4'1 = d^" ; (5.28) 



iW=0W-ivVi)-^aW; (5.29) 



(xi'^ll +2/3(1)) =0; (5.30) 



xni)+<i)=0; (5.31) 

f In other words, even if one starts with purely scalar linear perturbations as initial conditions for 
the second order theory, vector and tensor modes are dynamically generated jio) . 
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(1)T (1)T 

Ap 2j As 



(5.32) 



The parameters a'^^\ and dp"* of the gauge transformation can be obtained 



from (5.27), (|5.30D , and (5.31) respectively, while the transformed metric perturbations 
follow from ( js^ , (p.28[), (|5.29|), and ( 0^ ). 

Once these parameters are known, the transformation rules for the energy density 
/I or any other scalar, and those for the 4- velocity u^, follow trivially from ( 5.22 ), 
(5.24), and (5.25). 



5.2. Second order 

We now extend these well-known transformation rules of linear metric perturbations 
to the second order. 



5.2.1. General transformation Second order perturbations of the metric transform. 



according to (4.9), as 



(0) 



(5.33) 



This leads to the following transformations in the second order quantities appearing 
in (O)-(O): 



lapse perturbation 



2^(1) +2-^(1) 



"(1) 



shift perturbation 

4V'^^^a 



(1)' _L 

a 
(1) 



2a 



(1) (2^(1) 



^(1) 



"(2) + 7"(2) 



~(2) (2) 



„(1) 



(1)' 



a \ 
^(1) ( 

spatial metric, trace 



2(c.P + 2^c.«]-a«'+eP' 



0(2) = 0(2) + ^(1) 



2 0'(i) + 2-0(1) 



(5.34) 



(5.35) 



+ ^(1) i 2</'i'^ - - ^ (-40(1) + "(i)9o + + 4^a(i)) V^/3(i) 



i (2Kr) 



^1) 



+ ^(1) 



3 



(1) , t(l) I t(l)A rf":^ 
y ^ J S,'' i ''(I) 



(5.36) 



-"(2) - 3 V /3(2) ; 



19 



spatial metric, traceless part 



,~,(2) 



ii)' 



(1) 



,(1) Ck 



( 



2w 



(1) , 



- -(5„- 2cj 



fe' ^ 



(1) 



(1) 



(l)fe 



(1) , c(i) , efe.' 



(C-)+D.,/3(2)) . 



For the energy density /i, or any other scalar, we have from (4.9) 

From this we obtain, using (5.14): 

(5^/1 = S^^i + Ai(o)a(2) + a(i) (Ai(o)a(i) + M(o)"(i) + 2V 



€{!) (m(o)' 



(5.37) 



(5.38) 



(5.39) 



For the 4- velocity u'^, we have from (4.9): 



£ 



5(2) 



£i 



i(i) "(0) 



2£ea,<^«^ 



Using (5.8) and (5.15) this gives: 

^°2) = V{2) - ^"(2) - a(2) + "(1) 

- "(1) + ^(1) (2^('i) 



2 I 7,"' _ —7,0 

^ I «(1) ^^(1) 



a 1 1 / 
— a\- - a;- 
a ■' 



t'2 



+ ( 



"(2) = '"{2) - P(2) - «(2) + "(1) 



2 ( 



^(1) 



a 
a 



^(1) 



HI) ~ "'"Ml) 



+ 



(5.40) 



(5.41) 



^(1) (2«(i). - ^('d.) - ^(1). (2-fi) - eS)) + ^(i) (2^(1) + «'(!)) ; (5-42) 



for the time and the space components respectively. Again, the 4- velocity u'^ is subject 
to u^^u^J_ = —1, which gives (5.10); therefore (5.41) reduces to (5.34). 

5.2.2. Transforming from the synchronous to the Poisson gauge For this example, 
let us consider the simplified case in which only scalar first order perturbations are 
present as initial conditions for the second order problem. In the first order analysis 



presented above, this corresponds to having Xs ij = Xs ij 



(1)T _ 



7(1) 



,(1) 



X 



(1)J 

p ij 



^{i}± ~ ^' ^ii*^ ii^is 



0. The second order vector and tensor perturbations are 



however non- vanishing as the dynamical coupling of the modes makes them grow when 
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non-linear terms are considered in the evolution equation. We consider this restriction 
just for the sake of simplicity and because it describes a physically interesting situation. 
The more general transformation expressions follow straightforwardly from (5.34)- 
(5.37), (5.39), and (5.42). 

(2) 

Transforming from the synchronous to the Poisson gauge, the expression for Tpp 
can be easily obtained from (5.34), using ( ^.27 ) and the condition dl^-* = to express 
all the first order quantities in terms of P^^: 



(1) 



+ /^fi) I^P^P" + ^P^P'^ + 2P'(,] + + . (5.43) 



For Lu^^\ and (f)^^ we get: 



,(2) _ 



(5.44) 



— — ) ^(1) - -^(1) 



a 



For Xp \j obtain: 



(5.45) 



(5.46) 



Given the metric perturbations in the synchronous gauge, these constitute a set 
of coupled equations for the second order parameters of the transformation, a'^'^\ 
and dp\ and the second order metric perturbations in the Poisson gauge, 
(/)p , and Xp a- This system can be solved in the following way. Since in the Poisson 



P) - 



gauge d^Xp'ij = 0; we can use the fact that d'd^x'^\j = and the property d^d]'^' = 0, 
together with (5.46), to obtain an expression for from which /J^^^ can be 

computed: 



7(1) 



3^,(2):»i 
4 



72 -(l)r72. 



,(1).*^2«(1) 



.,f V72«(l)' 



1. 



V72v72fl(l) 



(5.47) 
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Then, using the condition d^Xp\j = ^^'^ substituting 13^"^^ we obtain an equation for 



+ + y/3S/3W - ^V^^(i)V^/3« (5.48) 

4 „, _o„nv 4 



FinaUy, using d^cu^^l = and substituting we get an equation for a^^^: 

- 2 (^24^) + - ^V2/3(i)) v2/3[i) - 2/?(^)/3«' . (5.49) 

Having obtained, at least imphcitly, all the parameters of the gauge transformation 
to second order, one can in principle compute the metric perturbations in the Poisson 
gauge from (5.43)-(5.46). 

Similarly, once the parameters are known, the perturbations in any scalar and 
4-vector, and in particular those in the energy density and in the 4- velocity of matter, 
follow trivially from (5.39)-(5.42). 

6. Conclusions 

In this paper we have studied the problem of gauge dependence in relativistic 
perturbation theory, considering perturbations of arbitrary order in a geometrical 
perspective. In fact, the problem itself is of a purely geometrical nature, dealing with 
the arbitrariness in the mapping between the physical spacetime and the background 
unperturbed one. Since no dynamics is involved, the formalism developed here can 
actually find application not only in general relativity, but in any spacetime theory. In 
considering a specific example, we have assumed a flat Robertson- Walker background, 
and derived the second order transformation between the well-known synchronous 
gauge , and the Poisson (generalized longitudinal) gauge discussed in . 

In linearized perturbation theory a gauge transformation is generated by an 
arbitrary vector field ^(i) , defined on the background spacetime, and associated with 
a one-parameter group of diffeomorphisms (a flow) : the gauge transformation of the 
perturbation ST of a tensor field T is then given by the Lie derivative £^^^^Tq of the 
background field Tq. However, in considering a gauge transformation from an exact 
point of view, we have found that it is not represented by a flow, but rather by a more 
general one-parameter family of diffeomorphisms. The question then was, how can we 
approximate the latter to a given order n? To this end, we have developed in section 
2 the necessary mathematical formalism. First, we have introduced certain families 
of mappings, dubbed knight diffeomorphisms of rank n, defined by ( ^.10 ). Then, in 



Theorem 2, we have proved that any one-parameter family of diffeomorphisms may 
always be approximated, to order n, by a knight diffeomorphism of rank n. This result 
(which confirms a claim in p[) is fundamental for gauge transformations of order n, 
as it guarantees that they are correctly represented by knight diffeomorphisms of the 
same rank. From the applicative point of view. Lemma 2 is thus all we need to have 
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a generating formula for the gauge transformation to an arbitrary order, (4.6). Since 
a knight diffeomorphism of rank n is basically the composition of n flows, and is thus 
generated by n vector fields ^(ij, . . . ,^(„), a gauge transformation of order n for the 
n-th order perturbation S"T of a tensor field T involves an appropriate combination 
of the Lie derivatives along ^(i), . . . , ^(„) of Tq, . . . , (5"~^T. 

Gauge transformations found their main application in considering the time 
evolution of perturbations of a given background spacetime: in a subsequent paper p5{ 
we shall look at second order perturbations of an Einstein de Sitter universe, comparing 
results in the synchronous and the Poisson gauges, thus applying the results presented 
in section 5. Beyond these applications, there are many topics that we have not 
touched upon here which, in a way or another, are related to gauge transformations, 
and become even more cumbersome in the non-linear case. We mention only a couple 
of them. We have implicitly assumed the applicability of the perturbative method; 
also, we have not considered the problem of eliminating spurious gauge modes. In 
particular, in our cosmological example, the synchronous gauge as defined in section 
5 should actually be regarded as a class of point identification maps [Q . 

A final issue we would like to mention is that of gauge-invariant quantities. In 
relation to this, we have first defined gauge invariance in an exact sense, and then given 
the conditions for the gauge invariance of a tensor field T to an arbitrary perturbative 
order n. However, we have not faced the problem of finding or constructing such 
quantities. In particular, in considering cosmological perturbations, it would be 
useful to have at hand a set of second order gauge-invariant variables defined d la 
Bardeen However, given the gauge-dependent metric perturbations and their 
transformation rules presented in section 5, the construction of such variables seems 
impractical. Moreover, it is far from obvious that a complete set giving a full second 
order description exists at all, as is the case at first order . Another possibility 

is to look for covariant quantities ^ : in order to be second order gauge- invariant 
these should vanish in the background and at first order. Assuming purely scalar first 
order perturbations, two examples are the magnetic part of the Weyl tensor and the 
vorticity of the 4- velocity of matter. Other second order gauge- invariant quantities 
can be defined by taking products of first order gauge-invariant tensors that vanish in 
the background. In particular, this is the case for scalars such as E^^E'-'-'^ , where E 
is the electric part of the Weyl tensor. Once again, it seems difficult that a complete 
set of such variables could even exist. Nevertheless, it is worth pointing out that 
quantities which are quadratic in first order gauge-invariant variables are useful in 
specific problems; for example they may intervene in the construction of effective 
energy momentum tensors of perturbations, which are important for the study of back 
reaction problems [p^ Another possible application of the formalism developed 
here can be the study of the gauge dependence of these quantities. 
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Figure 1. The action of a knight diffeomorphism generated by ^(i) and 4(2)- 
Solid Unes: integral curves of ^(i)- Dashed Unes: integral curves of C(2)- The 
parameter lapse between p and cf>^^^ (p) is A, and that from cj)^^^ (p) to 4>^^2^2 ° "^A^^ 
is AV2. 
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Figure 2. The action of a knight diffeomorphism of rank two, represented in a chart 
to second order. 




Figure 3. The action of a gauge transformation f>Xi represented on the background 
spacetime Mo by its second order approximation, generated by the two vector fields 
^(1) and ^(2) . 



4>„ o4>, (p) 

Figure 4. If X and Y do not commute, $ is not a flow on A4o- 



